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Upcoming fixed-target experiments designed to search for new sub-GeV forces will also have sen- 
sitivity to the never before observed True Muonium atom, a bound state of We describe the 
production and decay characteristics of True Muonium relevant to these experiments. Importantly, 
we find that secondary production mechanisms dominate over primary production for the long-lived 
2S and 2P states, leading to total yields an order of magnitude larger than naive estimates previ- 
ously suggested. We present yield estimates for True Muonium as a function of energy fraction and 
decay length, useful for guiding future experimental studies. Discovery and measurement prospects 
appear very favorable. 



I. INTRODUCTION 

A new generation of fixed-target experiments at Jef- 
ferson Laboratory and Mainz [DH] designed to search for 
new sub-GeV scale forces will also have unprecedented 
sensitivity to rare QED processes. Experiments such as 
the Heavy Photon Search (HPS) [3 will have the capacity 
to precisely identify displaced vertices in the e + e _ final 
state arising downstream of an electron beam scattering 
off a high-Z target. This opens the door to discovering 
True Muonium (TM), the QED bound state of a fi + n~ 
pair, and studying it for the first time. The decay and 
spectral characteristics of TM can then be used to fur- 
ther test properties of the muon, and used to study bound 
state physics, with non-perturbative analogues in QCD, 
in a calculable regime [3] ■ Studying TM is additionally 
motivated in light of long-standing discrepancies between 
theory and observations of (g — 2)^, as well as more re- 
cent discrepancies in the measured proton charge radius 
in muonic hydrogen [5] . 

In this short note, we compute total yields for TM 
production in fixed-target experiments like HPS, prop- 
erly including primary production, break-up, and exci- 
tation reactions of relativistic TM. HPS uses a forward 
peaked electromagnetic calorimeter to trigger on coin- 
cident charged particles or energetic photons produced 
in e~ -nucleus collisions. Then, using a relatively small 
silicon tracker, charged particles and any associated dis- 
placed vertices are identified within ~ 10 cm downstream 
from the target. This allows HPS to look for the pro- 
duction of rare QED states with long lifetime, such as 
TM. Our discussion will focus on the total production of 
triplet states l 3 Si, 2 3 Si, and 2 3 P 2 because these states 
eventually decay to e + e~ and can be detected by a silicon 
tracker. We do not aim to provide a precise prediction 
for the total TM production rates, but rather include all 
dominant effects of TM production and decay to help 
guide upcoming efforts to discover TM. In addition to 
the well-known primary production of TM, we show that 
secondary production via l 3 ^! excitations is the main 
source of 2 3 5i and 2 3 P 2 production. The 2S and 2P 
state are especially long-lived, so our finding is partic- 
ularly important for fixed-target experiments that can 



identify vertices of e + e _ pairs in the range of ~ 1 cm 
to several cm. 2 3 Si and 2 3 Pi produced via secondary 
mechanisms should be readily discoverable in an HPS- 
style experiment. 

In Section |Hj we describe our methods, and summa- 
rize our overall calculation and results for the yield of TM 
events relevant for upcoming fixed-target experiments. 
For illustration, and for ease of comparison with exist- 
ing literature, we present specific yield results for a Lead 
target. The results for Lead are also comparable to the 
results for other common high-Z targets like Tungsten or 
Tantalum. In Section |III[ we describe primary produc- 
tion calculations, followed by secondary mechanisms in 
Section |IV| We end with a short discussion regarding 
discovery and measurement prospects. 

II. PRODUCTION AND DECAY YIELDS 

For an HPS-style setup, the triplet l 3 5i, 2 3 Si, and 
2 3 P 2 will dominate the TM signal. The singlet config- 
urations decay to 77, which is very difficult to separate 
from QED backgrounds. Primary production through 
single- and 3- photon reactions [HJ [TO] dominate the total 
production of the l 3 5i state. 2 3 5i and 2 3 P 2 primary pro- 
duction is down by an order or magnitude compared to 
l 3 ^! 0HD]. Additionally, the dissociation for 2 3 S' 1 and 
2 3 P 2 is an order of magnitude larger than l 3 Si [TT1 112) . 
which makes primary production (followed by target es- 
cape) yields negligible. The leading source of 2 3 Si and 
2 3 P 2 production is through secondary reactions initiated 
by l 3 Si scattering into 2 3 P 2 , and 2 3 P 2 scattering or de- 
caying into 2 3 S\. In what follows, we will use primary 
production rates for l 3 Si, 2 3 S' 1 , and 2 3 P 2 previously cal- 
culated 0HU], along with l 3 5i 2 3 P 2 and 2 3 P 2 ->■ 2 3 5i 
transition and IS, 2S, 2P — > X excitation into all possible 
final states cross sections computed using the formalism 

of nuns]. 

We compute total TM yields in terms of the inci- 
dent number of beam electrons N e , as a function of dis- 
tance traveled through the target. In practice, differen- 
tial yields are a non-trivial function of x — E/Ei, earn . A 
convenient dimensionless unit of target thickness is the 
dissociation length (or more precisely, dissociation into 



2 



all final states length), l^s^x of the l 3 5i state. Thus 
we convert length I into z — -. — . Let N e (z) be the 

number of beam electrons (with energy Eb eam ) as a func- 
tion of distance z, and Ni$ t 2S,2P (dropping the spin la- 
bels) the average number of TM states. 

For the l 3 5i state, the yield as a function of distance 
is controlled by primary production and dissociation. To 
a good approximation, 
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where a(eT — > IS") is the cross-section for an electron to 
scatter off the nuclear target (Lead for illustration in the 
remainder of this paper) and produce the l 3 5i state, and 
a(lS —> X) is the cross-section for for the l 3 5i state to 
scatter and dissociate. For the 2 3 5i and 2 3 P 2 states, sec- 
ondary production mechanisms involving l 3 Si -J- 2 3 P 2 hz Sl ^ x is the mean free path for l 3 Si breakup. Blue is 
and 2 3 Si f-> 2 3 P 2 reactions are important, leading to a 
yield evolution well approximated by, 



FIG. 1: Relative abundance Yi = Ni/N e for 

i =l 3 S , i,2 3 S , i,2 3 P 2 TM states as a function of the distance 
traversed through target, in units of z = l/li^St^-Xi where 

l 3 5i, red is 2 3 Si, and dotted gold is 2 3 P 2 . 
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o(eT — > 25) and a(e~ — ¥ 2P) are the primary produc- 
tion cross sections for an electron to scatter of the nuclear 
target and produce a 2 3 5i and 2 3 P 2 states, respectively, 
while <t(25 — > X) and er(2P — > X) are the dissociation 
cross sections. er(25 — > 2P) is the cross-section for the 
2 3 5i state to scatter off the nuclear target into 2 3 P 2 , 
while er(2P — > 25) is the reverse reaction. These are just 
Boltzmann equations, and the yields quickly asymptote 
to constants as z exceeds the dissociation length set by, 
hs->X = A/N A pa(lS -> X). We will only be inter- 
ested in very thin targets, so beam spreading effects are 
ignored, and N e can be treated as constant. We have 
only included the states (and related transitions) that 
are populated by an amount larger than 5% of the to- 
tal TM yield. Triplet-singlet transitions are significantly 
suppressed [HJH2], so ^ e aDove evolution equations refer 
only to triplet states. Finally, we note that we actually 
solve this system as a function of the energy fraction x, 
which is obtained by replacing all cross-sections and yield 
factors by their differential form. We express final results 
in terms of relative yields, ^ = -h- as a function of 

x as well as total yields J* dx^-. 

Using the results of Sections [TTT] and |IV| the to- 
tal relative yields are shown in Figure [T] for Lead using 
a primary beam energy Eb eam — 6.6 GeV. For Lead, 
h 3 s 1 ^x ~ 4.4 x 10 -4 cm, or roughly 10 -3 radiation 
lengths. After a few dissociation lengths, the relative 



yields are constant, so in practice only the first 10~ 3 ra- 
diation lengths of the target are important for TM pro- 
duction. The asymptotic abundance of l 3 5i is controlled 
by the total production cross section and the dissocia- 
tion cross section. For 2 3 5i and 2 3 P 2 , the asymptotic 
abundance is controlled by a balance of l 3 5i — > 2 3 P 2 , 
2 3 5i -> 2 3 P 2 , and 2 3 5i,2 3 P 2 -> X dissociation reac- 
tions. In fact, the primary production terms for 2 3 5i 
and 2 3 P 2 are negligible compared to the secondary pro- 
duction terms, so in practice we can ignore them. After a 
few dissociation lengths, the relative abundances become 
approximately constant. For other high-Z materials, the 
relative abundance results are very similar - this is due 
to the approximate cancellation of the Z-scaling between 
the production cross-section and the dissociation reac- 
tions. 

In Figure [2j we show the relative yields as a function of 
x once three dissociation lengths li^Si-ix of target mate- 
rial have been traversed. The dissociation and excitation 
reactions are approximately independent of the TM en- 
ergy far away from threshold (see Section III). Thus, the 
x— dependence is controlled by the primary production 

terms 



rf<r(e~-»lS,2,S,2P) 
dx 



The peak at low-x is dominated 
by the 37 reaction of [TO], while the mild high-x peak is 
dominated by the single 7 brehmstrahlung-like reaction 

The primary trigger for HPS-style fixed-target experi- 
ments involves a cut on the energy fraction of the ob- 
served e + e~ pair, as this is required to remove the 
overwhelming rate of QED e + e~~ background at lower 
x = E e + e - 1 Ef, eam . An additional cut on the displaced 
vertex (typically 1-2 cm for 6.6 GeV beam energies) is 
required to remove all background. To present yield re- 
sults in a useful manner, we calculate the total yield of 
TM states as a function of cut, x c , and vertex cut, 1 C1 
assuming a beam energy of 6.6 GeV for reference. This 
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FIG. 2: Relative abundance Yi = Ni/N e for 

i =1 3 Si ,2 3 Si ,2 s Pa TM states at the steady-state configura- 
tion, as a function of energy fraction x = E / E beam . z = 10 
was chosen. Blue is l 3 Si, red is 2 3 Si, and dotted gold is 2 3 P2- 



is just, 



Y(x > x c , I > l c ) = I dx^2 



dYjjx) 
dx 



(4) 



where Ti is the lifetime of the i TM state. In Figure 
|3j we plot Y(x > x c , I > l c ) versus x c and l c for a Lead 
target and a beam energy of 6.6 GeV. These results are 
applicable so long as the target is thicker than ~ l X 3 
We note that the typical opening angle of the e + e~ pair 

is 5.8° x 



/ 6.6 GeV \ 
\ E beam )■ 



impact of adjusting beam energy is on the primary pro- 
duction cross-section. The total TM primary production 
cross-section as a function of x c is shown in Figure [4] 
While the overall growth of the cross-section with energy 
is mild (seen by looking at low x c ), the cross-section at 
x c > 0.8 sharply increases with energy. For example, the 
cross-section at x c = 0.8 increases by a factor of ~ 3 going 
from E beam = 2 GeV to E beam = 12 GeV. This is mainly 
due to the energy dependence of the bremstrahlung pro- 
cess, which both grows and more sharply peaks at high 
x as energy is increased [13] . Figure [5] illustrates the 
change in x-dependence with changing beam energy. 
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FIG. 4: Total primary production cross-section for all nS 
states as a function of energy fraction cut, x cu t- The target 
is Lead, and units are in pb. Blue (bottom) is for Ebeam = 2 
GeV, purple (middle) is for E bearn = 6 GeV, and gold (top) 
is for E beam = 12 GeV. 
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FIG. 3: Total yield of e + e events coming from the decays 
of l 3 5i,2 3 S , i, and 2 3 P 2 as a function of energy fraction cut 
(y-axis) and displaced vertex cut (x-axis in cm). The yield is 
given in units of 10 -18 -/V e , where N e is the number of elec- 
trons on target, and we've assumed that the target exceeds a 
thickness of a few l 1 3 Sl ^ x ■ 

In addition to boosting the decay displacements and 
altering the typical opening angle, the most significant 
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FIG. 5: Differential cross-section for primary production of 
all nS states as a function of energy fraction x. The target is 
Lead, and units are in pb. Blue (bottom) is for Eb ea m = 2 
GeV, purple (middle) is for E beam = 6 GeV, and gold (top) 
is for E beam = 12 GeV. 



III. PRIMARY PRODUCTION 

In this section, we summarize the essential primary 
production results used in this note. The differential pro- 
duction cross section for the Bremsstrahlug process from 
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9 for the n-th energy state is, 
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where x — -gp-, and Ef, eam is the beam energy. This re- 
action has a high-x logarithmic divergence regulated by 
m e , but in practice the Weizsacker Williams approxima- 
tion breaks down near iwl - m e /Ei )eam [SJ[T3]. Thus, 
we integrate up to the point where ^ — for computing 
the total cross-section. The variation of the cross-section 
integral that results by varying the cutoff used in this 
procedure is less than 5%. It is worth noting that the 
formulas for the total cross-section presented in [§] are 
not correct, and in fact run negative for beam energies 
of ~ 2 GeV. This appears to be the result of integrating 
the Weizsacker Williams form all the way up to x = 1. 
As a rough check on our results, we compare to a Mad- 
graph calculation for production of a heavy photon A' 
[15] . since they both have the same production mecha- 
nism. The comparison is based on matching the kinetic- 
mixing parameter e = g'/e (where g' is the coupling of 
A' to electrons), by requiring that the lifetime of a heavy 
photon of mass 210 MeV be that of TM. We get a 30% 
smaller cross-sections compared to our results, but this 
is within the systematic uncertainties of the comparison. 
For the 37 process, we use equation (10) from [TO] : 

1 „ a 7 ( ZaA\ 2 AB 
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MeV. We can readily see 
dependence compared to 



where B = 0.85 and A = 
that there is an additional 

the bremsstrahlung process, which yields proportionally 
larger primary production rates in heavier (high Z) tar- 
gets. For the total cross-section, we take the sum of the 
37 and bremsstrahlung reactions, ignoring interference. 
These reactions are peaked in different kinematic regions 
— Bremsstrahlung at high-x and 37 at low- a; — so we ex- 
pect that this approximation is good to 20 — 30%. In all, 
we expect that our primary production cross-sections are 
good to ~ 40%, based on uncertainties from atomic form 
factors, use of the Weizsacker Williams approximation, 
and interference effects. 



IV. SECONDARY PRODUCTION 

To compute secondary production effects (interaction 
of TM with matter), we use the formalism of [11], [12], 
and [2] . In this formalism, TM is treated to be initially 
in the rest frame, and matter is treated as relativistic 



structureless particles. This allows one to tractably con- 
sider excitations of TM, from which we can also calcu- 
late break-up cross-section. The main contribution to 
the cross-sections is from electric interactions, while spin- 
spin, para-ortho and magnetic (internal motion of the 
atom) interactions are smaller by a few orders of mag- 
nitude, and thus can be neglected. The more complete 
calculation of many of these effects is presented in |14) . 
but for atoms like Pb neglecting these leads to ~ 1% 
error. While the IS transition results are included in 
earlier literature (see [TTjrTS]), the relevant 2S and 2P 
transitions are included here for the first time. 

Before we calculate cross-sections, one has to calculate 
the transfer matrix elements: 



00 />TT />27T 

/ / x 2 Sm(0)e lclxCos(e) x 



ip nlm (x,6»,0)>™ (m (x,6», 



where "0's are the hydrogen-like atomic wave-functions 
for TM, q is the magnitude of the 3-vector part of the 
momentum transfer and x = - with a the Bohr radius 

a 

for TM. nl labels the incoming TM state (i.e. l 3 Si, 
2 3 Si... etc). We now have the electric differential cross- 
sections for transition into other bound states or break-up 
given by, 



&att = Z 2 —(l-F nl0 ' n ' l '°(q) 

IT \ 

X \\A(q,Z)\ 2 qdq. 



(7) 



Here A is the photon propagator in Thomas-Fcrmi- 
Molier form: 



A(q,Z)=4^ 



i=l 



> 2 + (3f ' 



(8) 



where A = with 



61=6.0,63 = 1.2,63 = 0.3, 
ax = 0.10, a 2 = 0.55, a 3 = 0.35 

The electric differential cross-sections for transition from 
the nl state to a specific n'l' state is, 



nl,n'l f 



= U-(-l) 



Z-Z' 



,a 2 1 



7r a* 



q\A(q,Z) 



x F 



nZ0,n'Z'0 (J) \2 Aq 



(9) 



To get total cross-sections, we integrate ^ and (|9| for 
q £ [0, 00). We can then calculate the dissociation cross- 
section by subtracting from the total "excitation cross- 
section" (integrated ^) the cross-sections for excitations 
into bound states (integrated Selected results are 

summarized in Table 1. These results agree with those 
listed in [TT]. All results are calculated using the same 
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formalism. The estimation of uncertainty for these exci- 
tation reactions is not straightforward, but a rough esti- 
mate can be derived by comparing to positronium results 
computed in the same formalism versus other methods. 
For example, the difference between results from [T5] and 
ours is below 10%. 



V. DISCOVERY PROSPECTS AND 
CONCLUSION 

We're now ready to assess the feasibility of TM dis- 
covery and possible measurements using upcoming fixed- 
target experiments. Typical experimental configura- 
tions consist of beam currents of a few hundred nA of 
continuous-wave electrons (in the case of HPS) with run 
times of O (months) and targets of thickness < 1% ra- 
diation lengths. Using our new calculation of total TM 
relative yields, we can express total TM production in 
terms of the number of electrons on target, assuming a 
target thicker than a few dissociation lengths of I 3 Si, 
which is typically < 10 -3 radiation lengths. For refer- 
ence values of / = 450 nA and 1 month of beam time, 
the total number N(x c , l c ) of TM produced with x > x c 
and I > l c is, 

AW c )^ c ,U(7.27xlO-) ( 45Q ///^J 

(10) 

For Lead (Tungsten is similar), we obtain total yields in 
the range of 200 to 600 events, depending on x c and l c . 
With total acceptance and reconstruction efficiencies in 



the range of 10 — 20%, this gives sizable detected yields. 

In this work, we have established that initial discovery 
of TM with the above experimental parameters should 
be possible. Moreover, TM yields should be large enough 
that measuring basic TM properties may be possible. In 
an upcoming paper [16] , a detailed analysis of target ma- 
terial, target configuration, and detector layout will be 
presented in order to define an optimal strategy for mea- 
suring TM production cross-sections, lifetimes, and dis- 
sociation rates. Additional work related to TM can be 
found in [T7H23]. 
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TABLE I: Cross-sections for excitations and dissociation (in 
cm 2 ) for Lead, "to e + e~" refers to dissociation of TM, while 
"X" refers to all final states inclusive. 
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